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Abstract. We analyze publicly available Fermi-LAT high-energy gamma-ray data and con- 
firm the existence of clear spectral feature peaked at E^ = 130 GeV. Scanning over the 
Galaxy we identify several disconnected regions where the observed excess originates from. 
Our best optimized fit is obtained for the central region of Galaxy with a clear peak at 
130 GeV with local statistical significance 4.5o\ The observed excess is not correlated with 
Fermi bubbles. We compute the photon spectra induced by dark matter annihilations into 
two and four standard model particles, the latter via two light intermediate states, and fit 
the spectra with data. Since our fits indicate sharper and higher signal peak than in the 
previous works, data favors dark matter direct two-body annihilation channels into photons 
or other channels giving only line-like spectra. If Einasto halo profile correctly predicts the 
central cusp of Galaxy, dark matter annihilation cross-section to two photons is of order ten 
percent of the standard thermal freeze-out cross-section. The large dark matter two-body 
annihilation cross-section to photons may signal a new resonance that should be searched for 
at the CERN LHC experiments. 
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1 Introduction 



If the existing cosmological dark matter (DM) [1] is a thermal relic consisting of weakly 
interacting massive particles, DM annihilations into the standard model (SM) particles should 
provide indirect evidence of DM in cosmic ray experiments [2]. In this scenario the first 
emerging signal of DM annihilations is expected to appear from Galactic regions with the 
highest DM density such as the centre of Galaxy or the largest DM sub-haloes. Very recently 
it was claimed [3, 4] that there is a 4.6<r (3.3cr) local (global) evidence of a monochromatic 
gamma-ray line [5, 6] with an energy ~ 130 GeV present in the Fermi Large Area Telescope 
(LAT) [7] publicly available data. The signal originates from the centre of Galaxy from a 
region obtained applying a signal-to-background optimization procedure on the gamma-ray 
data [4]. If the claim is true, this could be the very first evidence that the DM is of particle 
physics origin, representing a breakthrough both in cosmology and in particle physics. 

According to the analyses presented in ref. [4], fitting the Fermi-LAT excess with a 
narrow peak is more an assumption rather than a result. In fact, any sufficiently hard 
spectrum with sharp fall-off around 130 GeV, for example a box, would fit the data presented 
in figure 4 of [4] as well as the narrow peak. This would open a possibility to explain the 
excess with astrophysical sources, for example with an unknown mechanism associated with 
the Fermi bubbles [8]. As the best signal-to-background regions determined in [4] seem 
to overlap with the Fermi bubbles, such a qualitative connection is easy to come [9]. In 
addition, this would also allow to explain the observed excess with photon spectra from DM 
annihilations into standard model final states that produce significant amount of prompt 
photons in their decays [10]. Discriminating between those possibilities requires thorough 
study of the Fermi data. 

In this work we analyze the publicly available Fermi-LAT gamma-ray data in order to 
check and study independently the claim of E~ ~ 130 GeV excess. We analyze the Fermi- 
LAT 195 week ULTRACLEAN dataset using the kernel smoothing method for fitting that is 
independent of binning and is complementary to the sliding energy window method used by 
Weniger. To identify signal target regions we use data driven method similar to the one used 
by Fermi Collaboration in searches for DM sub- haloes [11]. We estimate the errors of our 
fits with the bootstrap method. Using Monte Carlo method we study what is the probability 
that the observed excesses are statistical fluctuations of background. We are interested in 
finding out whether the excess exists, what is its spatial distribution in the Galaxy, what is 
its spectrum, possible origin etc. 

Fitting the photon spectrum coming from the target regions identified by Weniger, we 
do confirm the existence a spectral feature in the analyzed data centered at Ej = 130 GeV. 
Knowing the energy of excess, we scan the data to find the regions where the excess comes 
from. We find that the excess originates from relatively small disconnected regions, the 
most important of them is the centre of Galaxy but several other regions exist. Away from 
the identified regions the excess disappears consistently with the expectations from DM 
annihilations. We fit the gamma-ray background at energies 20-300 GeV from data by cutting 
out the central signal region of the Galaxy. We obtain a perfect power-law fit with a power 2.6 
for the high-energy gamma-ray background. We then fit the spectrum from the central signal 
region and observe a clear peak at £L ~ 130 GeV with a local statistical significance 4.5a. 
Fits from other regions have local significances as high as 3.2<7. We comment on estimating 
the corresponding global statistical significances of those results. 

Our results indicate that the shape of the peak is even more pronounced than obtained 
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in [4]: we observe also a rise of the peak. This result disfavors any possible explanation to 
the observed excess with shallow spectra (including the box-like spectra) over the power-law 
background, and favors more peaked profiles that may be difficult to obtain from standard 
astrophysical sources. We study, classify and fit the possible model independent DM annihi- 
lation scenarios [12] into two and four standard model final states (the latter is assumed to 
occur via light intermediate states [13] that may also induce Sommerfeld enhancement of the 
annihilation cross-section). We find that, among those spectra, our results disfavor any other 
scenario but direct two-body annihilation into 77 or jX final states, where X is any mas- 
sive particle. However, other narrow peak-like spectra like from internal bremsstrahlung [3] 
or narrow boxes due to massive intermediate particles [10] are still allowed within present 
Fermi-LAT energy resolution. The DM annihilation cross-section to two photons should be 
of order ten percent of the standard thermal freeze-out cross-section. This result depends 
very sensitively on the unknown properties of the central cusp of our Galaxy and can there- 
fore be relaxed. Increasing the loop-suppressed DM annihilation cross-section to photons 
may require the existence of new resonances that should be searched for at the CERN LHC 
experiments. 

An important result of ref. [4] is showing that regions with optimal size, depending on the 
assumed DM halo profile, should be used for this type of search. Likely the reason why Fermi 
Collaboration did not observe the excess is that they looked at too large region [14]. Assuming 
DM halo profiles introduces theoretical bias into the analyses. In addition, we do not see a 
good reason why the low energy gamma-ray spectrum at energies 1 GeV < E~ < 20 GeV 
should be used to identify the best signal regions above > 100 GeV. To the contrary, 
we believe that using the low energy spectrum for these purposes may be misleading. For 
example, we observe a slight asymmetry in the low energy spectrum that explains the north- 
south asymmetry of the identified regions in [4]. Since for large regions the background is 
completely dominated by the Galactic disk, the requirement of good signal-to-background 
ratio cuts off the disk region. Thus the shape of the regions obtained in [4] must trivially 
be of a hourglass type - their overlap with the Fermi bubbles is most likely accidental. This 
coincidence lead some authors to speculate that the signal is associated with Fermi bubbles 
while the photons from the centre follow power- law background [9]. Our results show that 
the actual situation with the = 130 GeV excess is exactly opposite and the observed 
signal cannot be associated with the Fermi bubbles. 

Based on that criticism we use data driven method similar to [11] to search for signal 
regions, and we optimize their size ourselves. We find that the optimal signal region for the 
best statistical significance is covered by the radius of 3° in the Galactic centre. For larger 
regions the signal significance decreases due to larger background, for smaller regions, like the 
Reg5 in [4] , the signal decreases because of too small number signal photons from that small 
region. Thus we confirm the necessity of right choice of the signal regions. We also checked 
that the 3° region in the centre of Galaxy where the significance of the signal over background 
is maximized is consistent with the expectations from DM annihilations for Einasto profile. 
Our results show that the peak excess is concentrated to rather small regions. If the source 
of the 130 GeV photons is astrophysical, further investigation of those regions with different 
observation frequencies should reveal objects or processes that also produce the 130 GeV 
gamma-rays. However, if the origin of the excess is direct DM annihilation into photons, we 
may have identified the most dense DM regions of our Galaxy. This is anticipated result for 
the Galactic centre. However, for the other regions this claim must be confirmed by other 
experiments with more statistics because it is possible that we have observed just an upward 
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statistical fluctuation of the background [15]. 



2 Data analyses 
2.1 Data selection 

In the present analysis, we take into account 195 week (from 4 Aug 2008 to 18 April 2012) of 
data from the Fermi Large Area Telescope (LAT) with energies between 20 and 300 GeV. 1 We 
apply the zenith-angle cut 9 < 100° in order to avoid contamination with the earth albedo, 
as recommended by the Fermi LAT team. We also apply the recommended quality-filter cut 
DATA_QUAL= 1, LAT_CONFIG= 1, and ABS(ROCK_ANGLE)< 52. 2 We make use of the 
ULTRACLEAN events selection (Pass 7 Version 6), in order to minimize potential statistical 
errors. The selection of events as well as the calculation of exposure maps is performed using 
the 18 April 2012 version of ScienceTools v9r27pl. 3 For everything else we use our own 



In our analysis, we do not subtract the contribution from known point sources in the 
Fermi LAT data because the known astrophysical point sources are unlikely to be a problem 
at E > 80 GeV. However, we have checked that subtraction of known point-like sources do 
not affect our results in any way: the exclusion of known points sources do not affect the 
results when searching high-energy gamma-ray lines. 

2.2 Spectra estimation 

Calculating the spectra from observed events is practically a probability distribution esti- 
mation. Estimating probability distributions is a well-developed topics in statistics, and we 
can choose all the modern tools - kernel densities, adaptive kernels, smoothed bootstrap 
for point-wise confidence intervals. In [16] this method is applied to study the luminosity 
function of galaxies. In this paper we shall use the same method, briefly reviewed below. 

The simplest approach to find a observed spectra is the binned density histogram that 
depends both on the bin widths and the location of the bin edges. A better way to estimate 
a probability distribution is to use kernel smoothing (see, e.g. [17]), where the density is 
represented by a sum of kernels centered at the data points: 



The kernels K(x) are distributions (K(x) > 0, J K(x)dx = 1) of zero mean and of a typical 
with h. The width is an analogue of the bin width, but there are no bin edges to worry 
about. In the latter equation, we use the adaptive kernel estimation, where the kernel widths 
depend on the data, hi = h(Ei). The summation is taken over all data points (events). 

The kernel widths are known to depend on the density f(x) itself, with h ~ f(x)~ - 2 
(see, e.g. [18]): in regions of fewer data points, we use wider kernels. This choice requests a 
pilot estimate for the density that can be found using a wide constant width kernel. 

To estimate the spectra, we use B% box spline kernel: 



software. 




(2.1) 



B 3 (x) 



x - 2| 3 - 4|x - 1| 3 + 6\x 
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4|x + 1| 3 + \x + 2 



(2.2) 
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1 http:/ /fermi. gsfc.nasa.gov/ssc/data/access/ 

2 http: / /fermi. gsfc.nasa.gov/ssc/data/ analysis/LAT _caveats.html 

3 http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/ 
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This kernel is well suited for estimating densities - it is compact, differing from zero only in 
the interval x £ [—2, 2], and it conserves mass: Y2i Bs(x — i) = 1 for any x. This kernel is 
also close to Gaussian with a = 0.6. 

We used the logarithmic energy scale for our spectra estimation. For the pilot estimate, 
we used a wide kernel with the scale h = 0.1 (in logarithmic energy scale). This wide kernel 
leads very smooth distribution and is used to estimate the approximate event probability 
depending on energy. For the adaptive kernel widths, we adopted h = 0.03 (the typical 
uncertainties in the overall energy calibration of the Fermi LAT) as the minimal width (for 
the maximum pilot density) and rescaled it by the h ~ / p not (x) -0 ' 2 law. The event probability 
is higher in the low energies and drops at the high energies, leading to roughly twice as wide 
kernels for high energies in studied energy interval. 

If we choose the kernel width this way, we minimize the mean integrated standard error 
(MISE) of the density. We are also interested in the "error bars" , point- wise confidence limits 
(CL) for the density. This can be obtained by smoothed bootstrap [19, 20]. Here the data 
points for the bootstrap realizations are chosen, as usual, randomly from the observed data 
with replacement, but they have an additional smoothing component: 

E\ = Ej + hej, (2.3) 

where e is a random variable of the density K(x). 

We generated 10000 bootstrap realizations, using the adaptive kernel widths. We show 
the centered 95% confidence regions in our figures. 

Using the kernel smoothing method to estimate the probability distribution is appro- 
priate when the potential signal is weak or the number of data points is small. In this case 
the usual binning technique may hide the signal or introduce a false signal, depending on the 
bin width and bin location. The kernel smoothing method depends only on the number of 
data points and the used kernel widths. Whereby, the kernel widths take into account the 
uncertainties of the observed data points and this way the resulting distribution function is 
effectively unaffected of it. As a result, the estimated probability distribution reveals the 
true observed distribution as accurately as possible. 

3 Fits to data: signal versus background 

Our aim is to perform as model independent fits to the Fermi LAT data as possible. To 
achieve this goal we first find the high-energy gamma-ray background from data. For that 
we exclude the central Galactic region with a radius of 12° from data and fit all the data 
in the energy range 20 GeV < < 300 GeV. The choice to exclude the central region of 
Galaxy from the background fit is motivated by the expected signal in this region. However, 
actually the result is insensitive to doing that. The fitting procedure is described in the 
previous section. We obtain almost perfect power-law background estimate from data with 
the power 2.6. This result agrees with qualitative theoretical expectations and shows that the 
gamma-ray background in this energy region is induced by scattering of high-energy diffuse 
protons. We use the data-fitted background in all our computations. 

Logically the first step towards more thorough analyses is to check the claim of the 
existence of E^ = 130 GeV spectral feature in the data. For that we first choose the Reg3 
identified in ref. [4] as our signal region and fit the gamma-ray data from that region as 
described in the previous section. The result is shown in figure 1 where we plot the resulting 
gamma-ray spectrum (red solid line) as a function of photon energy for the high-energy region 
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Figure 1. Estimated high-energy gamma-ray spectrum originating from Reg3 of ref. [4] together 
with 95% CL error band as a function of photon energy. Data fitted background (black solid line) is 
also shown, the pawer-law spectrum with power 2.6 (dotted line) is plotted for comparison. The blue 
dashed lines show 95% CL for statistical fluctuations of the background. 



20 GeV < E~ < 300 GeV. The 95% CL error band is shown with a grey band around the best 
fit. The background obtained from data (black solid line) and the perfect power-law spectrum 
with power 2.6 (dotted line) are also shown. The total number of high-energy photons and 
the number of photons with energies 120 GeV < E 7 < 140 GeV coming from this signal 
region are presented in table 1. The excess at E~ = 130 GeV with statistical significance 4 
3.6a is clearly visible. We observe that the excess has a more pronounced peak-like shape 
than that presented in figure 4 of ref. [4]. However, the overall flux we obtain is in good 
agreement with the gamma-ray flux obtained in ref. [4]. 

The question one may ask is whether the observed excess in figure 1 could entirely be 
due to statistical fluctuations of the background. Another way to estimate the significance of 
the peak is to generate background fluctuations with Monte Carlo for the number of photons 
in the target region. To do that we generated 10000 Monte Carlo realizations where the 
photon energy distribution follows the observed background: the high-energy background 
shape is assumed to be the same over the sky. For every realization we calculated the spectra 
as described in previous section. Based on 10000 Monte Carlo realizations we extracted 95% 
confidence level limits for statistical fluctuations of the background. The results are presented 
in figure 1 with blue dashed lines. While all other small excesses and deficits around the 
background are consistent with expectations of fluctuations, the signal at 130 GeV clearly 
exceeds 95% CL for background. 

Having confirmed the previous claims, the next question to ask is from which region 
of the Galaxy the photons in the peak come from? For illustration of data we first plot 
the distribution of photons in the energy range 120 GeV < E 7 < 140 GeV, denoted with 
blue dots, in the left panel of figure 2. Thus the figure is actually a Fermi photograph of 

4 The statistical significance is estimated using the bootstrap realizations as described in section 2.2. The 
significance is given respect to the calculated background of the spectrum. 
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our Galaxy in this energy range. As expected, we observe that most of the photons in that 
energy range come either form the very centre of the Galaxy or from the Galactic disk area. 

In order to find the spatial origin of the 130 GeV excess we scan over the Galaxy 
choosing target regions with varying radii and fitting data coming from those regions as 
described before. Depending on the location we obtain either an excess or a deficit of the 
signal. At the same time we also compute 68% CL for the background fluctuation for that 
region as just described: the background shape is assumed to be constant over the sky, and 
the CL depends only on the number of observed photons. We define relative signal intensity as 
a ratio of the signal excess/deficit and the 68% confidence level background fluctuation limit 
(although it is in units of background fluctuation as, we do not want to confuse the defined 
relative intensity with the statistical significances of the signal calculated with bootstrap). 
Scanning over the Galaxy we obtain a distribution of the relative signal intensity of the 
excess/deficit in the Galaxy. To eliminate very large statistical fluctuations, we calculate the 
signal intensity only for regions where the number of photons in energy range 20-300 GeV 
is larger than 80. 

We plot in the right panel of figure 2 the resulting distribution of relative signal in- 
tensity as presented by the colour code. The pink background is due to regions with too 
low photon flux to obtain statistically meaningful results. As seen in the figure, the signal 
with highest significance originates from the centre of Galaxy. This region is centered at 
(I, b) = (— 1°,— 0.7°), called "Central" region in the following, and has a radius 3°, drawn 
with a white circle in figure 2. The total number of high-energy photons and the number of 
120 GeV < E<y < 140 GeV photons coming from this signal region is presented in table 1. 
However, there exist other regions, spatially well separated from the centre, that also exhibit 
large 130 GeV gamma-ray excess over the background. The most significant of them, with 
the same radius, is located at (I, b) = (—10°, 0°), called "West" region in the following, and is 
also shown in the figure. Some other possible signal regions are all listed in table 1. Presently 
statistically significant fits are obtained only for the first two regions, but with more Fermi 
statistics the other regions may become relevant too. 

One can see in figure 2 that the regions with excesses and the regions with deficit 
of the signal are not in balance - the excess dominates (considering the degree of dark- 
ness/brightness of the regions). The deficit almost never exceeds 2a level and is in good 
agreement with the expectations from statistical fluctuations of the background. At the 
same time, there exist regions in which the observed excess exceeds the 2a level. 

It is clear from figure 2 that the excess of photons with energy around 130 GeV does not 
originate from Fermi bubbles. Firstly, there is no spatial correlation between the signal excess 
and the Fermi bubbles. Secondly, whatever is the physical mechanism creating the 130 GeV 
excess, this mechanism must be at work in several regions of the Galaxy. If the origin of the 
excess is astrophysical, it should be possible to observe those astrophysical objects/processes 
in the identified regions with other methods. Any such a mechanism must also explain why 
the observed excess is a peak, that might be difficult in the case of standard astrophysical 
processes. If, however, the origin of the 130 GeV peak is DM annihilations, figure 2 shows the 
distribution of the most dense DM sub-structures in the central region of our Galaxy. Notice 
that the dark centre of the Galaxy does not exactly coincide with the galactic coordinate 
origin. 

The fits to high-energy gamma-ray data originating from the Central and West sig- 
nal regions are plotted in the left and right panels of figure 3, respectively, using the 
same notation as in figure 1. The Central region exhibits an excess with local statisti- 
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Figure 2. Left: a Fermi "photograph" of our Galaxy in gamma-rays with the energy 120 GeV < 
E 1 < 140 GeV. Fermi data is shown with blue dots. The region Reg3 of ref. [4] (dashed black curve) 
and the Fermi bubbles (red areas) are also shown for illustration. Right: distribution of relative signal 
intensity of 130 GeV photons in the Galaxy. The green circles denote the signal regions that provide 
the excess with highest statistical significance; grey circles denote other regions showed in table 1; 
green dot mark the assumed centre of the Galaxy. 

cal significance of 4.5a. To estimate the global significance we consider the trials factor 
~ (energy range under study) /(width of the line), which reduces the global significance to 
4.0a. The angular size of the region is placed by the null hypothesis of the annihilation 
signal from the DM halo profile. The difference of the signal in case of the Einasto and 
NFW profiles and also the off-axis location of the Central region can slightly reduce the 
global significance further. The fit to West region shows a clear peak at 130 GeV with local 
statistical significance of 3.2a. We have also fitted the signal from other bright regions in 
figure 2 that all show an excess peaked at the same photon energy, £L = 130 GeV. Those 
are listed in table 1. The trial factor of the those regions is larger due to large region of 
sky scanned for the signal, allowing for the possibility that the excess in those regions is an 
upward fluctuation [15]. This issue must be clarified by future experiments. 

Based on the model independent results presented in figures 1-3 and in table 1, we 
conclude that, whatever is the physics origin of the excess, its significance is high, it has a 
clear peak shape, and it comes from the small region in the Galactic centre and possibly from 
several other small regions we have identified. 

4 Fitting DM annihilation spectra 

4.1 Comparison of different annihilation channels 

It is appealing to attempt to explain that the observed gamma-ray excess with DM anni- 
hilations in the Galactic centre and in the most dense DM sub-haloes around the centre of 
our Galaxy. Our approach is model independent as described in refs. [2, 12]. We assume 
that the DM particles annihilate into two SM particles, DM + DM — > SM + SM, where 
SM = 7, e, /i, r, q, W, Z, h, where q denotes any light or heavy quark. The final state 
SM particles decay and/or radiate photons and light fermions from the final state radia- 
tion. In addition, we also allow DM annihilations into two light hypothetical finals states, 
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Table 1. Identified signal regions in the Galaxy, number of photons in the two energy intervals and 
the statistical significance of excess in those regions. The radii of regions are all 3° (except for Weniger 
Reg3). 



Region 


I (deg) 


b (deg) 


iV 7 (20-300) GeV 


iV 7 (120-140 GeV) 


significance 


Weniger Reg3 






3298 


65 


3.6(7 


Central 


-1 


-0.7 


818 


27 


4.5cr 


West 


-10 





726 


21 


3.2cr 


East 


17 


-3 


481 


14 


2.7cr 


North 


-7 


16.5 


109 


4 


1.6(7 




Figure 3. Best fits to high-energy gamma-ray data for the Central (left panel) and West (right panel) 
signal regions presented in table 1, together with 95% CL error band as functions of photon energy. 
Background fitted from data is also shown (black solid line), the power-law spectrum with power 2.6 
is plotted for comparison (dotted line) . The blue dashed lines show 95% CL for statistical fluctuations 
of the background. 



DM + DM — > V + V, that decay as V — > 77 , ee, \x\i. Those particles have been postulated 
to exist [13] in order to explain the DM annihilations to lepton and not quark final states. 
In addition, those light particles may induce Sommerfeld enhancement of the annihilation 
cross-section. Both those features are needed to explain the PAMELA anomaly of positron 
flux with DM annihilation scenarios. Assuming those annihilation channels we have com- 
puted the resulting decay /hadronization chains and final state radiation with PYTHIA [21], 
and obtained the resulting prompt photon spectra. We use those prompt spectra to fit the 
observed gamma-ray excess. We neglect photons from inverse Compton scatterings between 
possible charged annihilation products and the Galactic and CMB photons as this spectrum 
is too distributed to explain the observed peaked excess. 

To study which DM annihilation scenarios can explain the observed excess best we 
estimate the goodness of the fit as follows. To fit the peak we choose the signal energy range 
to be 50 GeV < E 1 < 200 GeV and divide it into n = 100 bins. We compute \ 2 according to 

" / f th fobs \ 2 

i 1 

where f° hs and (7° bs are the observed flux and its error from our fits to data for a bin i, and 
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Figure 4. Best fit to data from the Central region with DM annihilation spectra into 27 (left panel) 
and to 2V — > 47 (right panel), presented with dashed blue lines. The rest is as in figure 3. The latter 
case can fit data only if large statistical fluctuation of the spectrum occur due to the present limited 
statistics, as demonstrated in the right panel with green dotted lines. 



/* is a theory estimate of signal plus background computed for every annihilation channel. 
The theoretical spectrum is calculated the same way as observed spectrum (see section 2.2) 
to have comparable spectra. For every annihilation channel we find the minimal reduced 
X 2 /n that gives the best fit of that theoretical model to data. 

Contrary to our initial intuition based on ref. [4] results, almost all the resulting photon 
spectra from DM annihilations turn out to be too soft to explain the observed gamma-ray 
excess. Our best fit (for the Central signal region) for the annihilation channel DM + DM — > 
7+7, denoted by blue dotted line, is presented in the left panel of figure 4 together with the fit 
do data as in figure 3. The best fit x 2 /n = 0.7 is obtained for the DM mass Mdm = 130 GeV. 
It was shown in ref. [4] that DM annihilations into a monochromatic gamma-ray line can fit 
the data. Our results agree with this conclusion. 

The second hardest photon spectrum is obtained for DM +DM — > V+V — > 47 channel, 
all the remaining studied channels produce so broad photon spectrum that cannot fit the 
observed peak. We present the best fit for this channel with x 2 / n = 2.7 in the right panel 
of figure 4 for the DM mass -Mdm = 145 GeV. This annihilation channel can reproduce the 
fall-off of the peak but predicts flat box-like distribution for smaller than the peak energies. 
The reduced x 2 /n is 2.7 compared with 0.7 for the 27 channel and the fit is clearly worse. It 
is still possible that, due to limited statistics, we observe significant down-ward fluctuation 
of the flux at photon energies ~ 100 GeV. With the present statistics the probability this 
to happen is at the level of a few percent. Examples of such fluctuations are presented in 
the right panel of figure 4 with green dotted lines. If we observe the down-ward fluctuation, 
this should disappear with more data and the observed peak will be replaced by a box-like 
spectrum. This scenario looks, however, unlikely. 

In the previous computations we assumed almost massless intermediate particles V. 
However, the width of the box-like spectrum depends on the mass difference between the 
DM mass and the intermediate state mass [10]. For the fine tuned case Mdm ~ My the 
box-like spectrum shrinks into a narrow peak-like spectrum, and a good fit to data becomes 
possible. 
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4.2 Fitting annihilation cross-section 

Up to now we have made no assumption about the DM distribution in the main DM halo 
of the Galaxy. However, if we want to estimate the annihilation cross-sections that fit the 
signal, we must compute the absolute flux of gamma-rays coming from the signal regions. 
Here we work only with the Central signal region. We assume two different halo profiles, the 
Einasto profile [22-25], 

(4.2) 

with a = 0.17, p s = 0.079, and r s = 20 and the Navarro-Frenk- White (NFW) profile [14, 26], 

PNFWW = p.- ( 1 + -) , (4.3) 

with p s = 0.33 and r s = 20. The profiles have been normalized to the local DM density 
Pdm = 0.4 GeV cm -3 at the Solar system [27, 28]. The gamma-ray flux from Central signal 
region is calculated according to formalism presented in [2]. To estimate the annihilation 
cross-section, {av), we let it to vary and compute the corresponding x 2 /nof the flux according 
to (4.1). The minimal x 2 / n gives the best fit value for {av). To estimate the confidence limits 
for cross-section, we use the bootstrap sample calculated from the observed events. To every 
bootstrap realization, we fit the model spectra as described in section 4.1 and find the best- 
fitted cross-section, {av). All bootstrap realizations give the distribution for {av), where we 
extract the 95% confidence level limits given in table 2. 

The results for 2DM — > 27 and 2DM — > 2V — > 47 annihilation channels with very light 
V are presented in table 2. They correspond to the cases plotted in figure 4. The annihilation 
cross-sections to photons are required to be large, of order ten percent of standard thermal 
cross-section, with quite large errors. Within errors our result for annihilation cross section 
agrees with the result obtained by Weniger. 

The resulting DM annihilation cross-section is much larger than expected in most of 
DM models. The result can be explained by cross-section enhancement mechanisms like 
the Sommerfeld enhancement or resonances in the annihilation process. We will elaborate 
on the latter in the next section. The enhancement of the flux can also be explained by 
non-standard DM cusp in the centre of Galaxy. In this case the annihilation cross-section 
can be small with the price of making the central cusp DM density bigger that predicted by 
the profiles we assume. If the main halo has substructures in the central region as favored 
by our results the presented cross-sections can have larger uncertainties as presented here. 
The over-dense substructures can explain the enhancement of the observed flux. In this case, 
assuming some theoretical value {av) for the annihilation cross-section, our results allow to 
compute the required DM density in the cusp of Galaxy and in the most dense DM sub-haloes 
that induces the measured flux. 

5 Discussion of the results 

While we do confirm the existence of the 130 GeV excess in Fermi LAT data, our results differ 
from the ones presented by ref. [4] in some important aspects. Firstly, our search strategy for 
finding the most sensitive signal regions in the Galaxy gives significantly different results than 
obtained by Weniger. We do not use the low energy gamma-ray data, 1 GeV < E y < 20 GeV, 



pmn(r) = p s exp 
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Table 2. The annihilation cross-sections for Central region given in units of the standard thermal 
cross-section, (av) — 3 x 10~ 26 cm 3 s _1 together with 95% confidence level limits. 



Channel 


(<Tf ) for Einasto 


(<tu) for NFW 


2DM -> 2 7 
2DM -»• 2F -»• 4 7 


0.09 ±0.05 
0.13 ±0.09 


0.17 ±0.09 
0.25 ±0.17 



to determine the high-energy gamma-ray target regions. At such low energies the background 
may be distorted by astrophysical sources like Fermi bubbles and, indeed, we observe north- 
south asymmetry in the low energy data. This explains why the Weniger target regions are 
asymmetric. In addition, Weniger's signal regions similarity with Fermi bubbles is most likely 
accidental but misleading. At the end, the fit to data from those regions gives quite significant 
excess, see figure 1. Instead, we used a data driven method to find the best target regions and 
we fitted the background from high-energy data by excluding the centre of Galaxy. With this 
procedure we identified small regions where the signal significantly exceeds the background, 
presented in figure 2 and table 1. Similar method was used by Fermi Collaboration to search 
for DM sub-haloes. In the context of our approach, the chosen target regions may correspond 
to DM sub-haloes of our Galaxy. 

Secondly, our results show that the excess has clearly a shape of the peak while Weniger's 
data analyses may also be consistent with a box-like excess. This result narrows the pos- 
sible astrophysical explanations to the excess. This result also disfavors all DM annihi- 
lation modes but the ones to photons with narrow peak-like spectrum (including internal 
bremsstrahlung [29-31] spectrum and narrow box- like spectrum). At the moment the width 
of the allowed peak is mostly determined by the Fermi LAT energy resolution and quite poor 
statistics. 

We found that the DM annihilation cross-section (av) should be larger than naively 
expected for loop-induced annihilation processes. Indeed, typical suppression for one loop 
processes involving weakly interacting particles is of order 10 -4 compared to the typical tree 
level processes. There are two possible ways to explain this result. Firstly, the DM annihi- 
lation cross-section into 77 may be small, but the DM halo models we have used, Einasto 
and NFW, fail to describe the central cusp of our Galaxy. This is possible since N-body 
simulations cannot predict the central region precisely. Since the cross-section calculation 
depends on the halo properties, there is always large theoretical uncertainty related to that. 
Secondly, the annihilation cross-section into photons might be, indeed, enhanced. Thus the 
model building of DM theories should concentrate on 77 or 7X annihilation modes and on 
the enhancement of those annihilation cross-sections to the observed level. 

Generically one expects DM direct annihilation cross-section to photons to be orders 
of magnitude smaller than is needed to explain the observation. However, examples exist in 
which loop level [32-37] or anomaly induced [38, 39] cross-sections to photons are enhanced. 
Motivated by the recent LHC results that at least one light fundamental scalar, the Higgs 
boson, likely exist with a mass 125 GeV, and no low scale supersymmetry nor Kaluza-Klein 
particles exist below 1 TeV, one can construct most plausible scenarios of thermal relic DM 
that have large annihilation cross-sections to 27 channel. 

For example, one well motivated possibility is that the DM consists of dark scalars 
that couple to the SM fermion via extended Higgs sector with non- vanishing (or dominant) 
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coupling only to the top quark. In this case the correct DM thermal relic abundance is induced 
via s-channel annihilation process for DM masses somewhat below the top quark mass [36]. 
This is a generic prediction of the set-up not related to any model building detail but to the 
topology of freeze-out process and to the top mass. Today, when the DM is essentially at 
rest, the only kinematically allowed annihilation channel to the SM particles is opened by 
the top loop induced coupling to two photons (that dominates over the jZ, ZZ final states). 
The kinematically allowed annihilation channels via virtual top-quarks imply 6-particle final 
states and are negligible. This, as argued before, must be enhanced. The enhancement can 
occur due to an accidental resonance so that 2Mdm ~ Ms, where Ms is the mass of new 
unobserved scalar. In this case all the parameters of the scenario are essentially fixed. The 
DM mass 130 GeV is consistent with thermal freeze-out, and new particle with non-standard 
couplings is predicted to exist with a mass Ms ~ 260 GeV, depending somewhat on how 
close to the resonance the annihilation process must be. A drawback of the scenario is that 
the DM-nucleon scattering cross-section is predicted to be small, explaining the absence of 
signal in DM direct searches at XENON100 [40] . 

This scenario has already been realized in ref. [36] for the 7/1 final state, where h is the 
SM Higgs boson, in the context of Randall-Sundrum framework with extra Z' boson. For the 
final states jX the 130 GeV gamma-ray excess implies different DM mass that depends on 
the mass of X via E 1 = M DM (1 - M|/4Mg M ). Apart from the mass shift, to enhance the 
gamma-ray signal to the indicated level, all this type scenarios should include new resonance 
at 2Mdm that can be searched for at the LHC experiments. 

6 Conclusions 

We have analyzed Fermi LAT publicly available data collected during 195 weeks and confirm 
the existence of gamma-ray excess peaked at £L = 130 GeV. The excess originates from 
disconnected spatial regions, presented in table 1. The excess is not correlated with Fermi 
bubbles. The strongest signal over background comes from the central region of the Galaxy 
and has local statistical significance of 4.5cr (global 4. Oct after trials factor correction for 
energy), but the signal from the other regions shows also up to 3a excess (in local significance). 
According to our fits the excess is narrower and higher than shown before. Leaving aside the 
possibility that the excess is an instrumental artefact, our results show that the mechanism 
of generating such an excess must be at work in several regions of our Galaxy. It might be 
difficult to explain the sharp gamma-ray peak with standard astrophysical processes. It is 
more appealing to assume that the excess originates from DM annihilations. In this case our 
results imply that we have identified the most dense DM substructures of our Galaxy - the 
central cusp and possibly some DM sub-haloes, although the latter claim must be confirmed 
with better statistics in future experiments. 

Assuming the DM annihilation scenario, we have computed photon spectra from DM 
annihilations into two SM particles and into two light bosons V that decay to two photons 
or leptons using PYTHIA. We find that only the DM annihilations into 77 or 7X final 
states, where X is any massive particle, can fit data well. All other spectra, including 
the 47 one from light intermediate V states, provide significantly worse fits to data. The 
exceptions are the possibilities that the V mass is almost degenerate with the DM particle 
mass since then the box-like spectrum shrinks into a peak, and the internal bremsstrahlung 
like spectra. We obtain that the DM annihilation cross-section into photons must be of 
order ten percent of the standard thermal freeze-out cross-section, see table 2. However, 
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this result is correct only if the Einasto or NFW halo profiles predict the properties of the 
central cusp exactly. Conversely, assuming some theoretically motivated DM annihilation 
cross-section to photons, our results allow to calculate the DM density in the cusp of Galaxy 
and in the most dense DM sub-haloes. We have sketched a generic thermal relic DM scenario 
that, independently of model details, should produce the observed DM relic density and the 
enhanced DM annihilation cross-section into 77 today. This scenario predicts a new resonance 
close to 2Mdm with specific couplings that should be searched for at the LHC experiments. 
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We have updated the fits to Fermi-LAT publicly available gamma-ray data from the 
Galactic centre presented in Ref. [41] using 218 weeks data and new improved Fermi-LAT 
energy resolution [42]. The new result is presented in Fig. 5 that shows a fit to the gamma- 
ray flux as a function of energy together with 2a error band. Compared to Figs. 3. and 
4. of Ref. [41], the most important new feature is the presence of a double peak in the 
130 GeV excess. While the previous Fermi-LAT energy resolution did not allow us to see the 
fine structure of the excess, the improved one clearly indicates for a double peak structure, 
confirming similar claims made in Ref. [43]. 

The double peak at the same energies, 110 GeV and 130 GeV, is also observed in the 
gamma-ray excess from nearby galaxy clusters [44] , suggesting that the two signals originate 
from the same source. The presence of double peak is a generic prediction of Dark Matter 
annihilation pattern in gauge theories, corresponding to 77 and final states. Thus the 
two seemingly unrelated gamma-ray spectra, from the Galactic centre and from the galaxy 
clusters, favour the particle physics origin of the excess over any astrophysics origin. 

We finally note that the double peak is not present in the gamma-ray spectrum from 
Earth Limb [45, 46]. 




20 100 130 200 

E [GeV] 



Figure 5. Updated fit to Galactic centre data using 218 week data and improved Fermi-LAT energy 
resolution [42]. Vertical lines mark the 110 GeV and 130 GeV energies. 
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